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impedance.  Thus  an  equivalent  circuit  is  required  which  can  predict  the  effect  of  eddy  currents  on 
the  efficiency,  the  electrical  impedance,  the  transmitting  responses,  and  the  receiving  response  of  a 
m^netostrictive  transducer.  *n»ls>eport  presents  the  derivation  of  an  equivalent  circuit  of  a 


magnetostrtctive  transducer  In  the  presence  of  eddy  currents  using  an  exact  eddy-current  theory. 
However,  use  of  the  exact  theory  can  be  time  consuming.  Lumped  equivalent  circuits  are  also  pre- 
sented which  use  approximations  to  the  exact  theory  to  simply  and  accurately  predict  the  efficiency 
and  electrical  impedance  of  a transducer  with  eddy  currents  within  a frequency  range  of  practical 
Interest  near  the  clamped  characteristic  frequency. 
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THE  EQUIVALENT  CIRCUIT  IN  THE  MOBILITY  ANALOGY 
OF  A MAGNETOSTRICTIVE  TRANSDUCER 
IN  THE  PRESENCE  OF  EDDY  CURRENTS 


INTRODUCTION 

Recent  research  into  the  highly  magnetostrictive  rare  earth  iron  compounds  has  renewed 
an  interest  in  magnetostrictive  transducers.  A typical  magnetostrictive  transducer  is  designed 
with  the  aid  of  an  equivalent  circuit.  The  equivalent  circuit  is  an  electrical  representation  of 
the  differential  equations  that  govern  the  motion  of  the  transducer.  The  equivalent  circuit  is 
used  to  predict  the  transmitting  responses,  receiving  response,  electrical  impedance,  and 
effioieney  of  the  transducer.  .....  .... 


Ohmic  losses  generated  by  eddy  currents  are  a primary  loss  mechanism  in  metallic 
magnetostrictive  transducers.  These  ohmic  losses  reduce  the  efficiency  and  alter  the  electri- 
cal impedance  of  a magnetostrictive  transducer.  Thus  an  equivalent  circuit  is  needed  which 
can  simply  and  accurately  predict  the  effect  of  eddy  currents  on  the  efficiency  and  the  elec- 
trical impedance  of  a magnetostrictive  transducer. 

The  active  elements  of  a magnetostrictive  transducer  are  typically  in  the  shape  of  a 
thin  rod  or  a stack  of  thin  laminations.  For  these  standard  configurations  it  would  be  highly 
desirable  to  be  able  to  predict  the  efficiency,  receiving  and  transmitting  responses,  and  elec- 
trical impedance  as  a function  of  eddy  current  loss,  which  is  controlled  by  the  rod  diameter 
or  lamination  thickness. 

The  first  purpose  of  this  report  is  to  show  how  the  equivalent  circuit  of  a magnetostric- 
tive transducer  in  the  presence  of  eddy  currents  is  derived.  Since  the  equivalent  circuit  based 
on  the  exact  theory  can  be  time  consuming  to  use,  the  second  purpose  of  this  report  is  to 
present  useful  approximations  to  an  exact  eddy-current  theory  that  can  simply  and  accu- 
rately predict  the  efficienty  and  electrical  impedance  of  a magnetostrictive  transducer  in  the 
presence  of  eddy  currents.  The  third  purpose  is  to  discuss  the  effects  of  eddy  currents  on  the 
receiving  and  transmitting  responses  of  a magnetostrictive  transducer.  The  fourth  purpose  is 
to  discuss  the  effect  of  magnetomechanical  coupling  on  eddy  current  loss. 


THEORY  AND  DISCUSSION 

Eddy  Currents  in  a Rod  and  a Thin  Sheet 

When  an  alternating  current  is  applied  to  a solenoid,  a circumferential  electric  field  is 
generated  by  the  time-changing  flux.  If  a conducting  rod  is  placed  in  the  solenoid,  the  elec- 
tric field  generates  circumferential  currents  in  the  rod  (Fig.  1).  These  cun-ents  are  in  a direc- 
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Fig.  1 — Eddy  current*  in  • conducting  rod.  The 
current  i(f)  i*  increasing  with  time. 


tion  so  as  to  generate  a magnetic  field  that  opposes  the  applied  magnetic  field.  These  circu- 
lating currents  are  known  as  eddy  currents.  In  a rod  the  magnitude  of  the  eddy  currents 
increase  in  direct  proportion  to  the  radial  distance  from  the  center  of  the  rod.  Thus  near  the 
center  of  the  rod  the  opposing  magnetic  field  generated  by  the  eddy  currents  is  the  largest, 
since  all  the  eddy  currents  encircle  this  axis.  This  effectively  shields  the  center  of  the  rod 
from  the  applied  magnetic  field,  leading  to  the  well  known  skin  effect.  Since  the  rod  in  the 
solenoid  has  a finite  resistivity,  the  eddy  currents  also  generate  an  ohmic  loss.  The  shielding 
effect  and  ohmic  loss  generated  by  the  circulating  eddy  currents  alter  the  electrical  imped- 
ance of  a solenoid  with  a metallic  core.  This  subsection  will  present  the  derivation  of  the 
electrical  or  blocked  impedance  of  a solenoid  with  a core  composed  of  a rod  or  a thin 
sheet  in  the  presence  of  eddy  currents. 

The  penetration  of  a magnetic  field  into  a conductor  is  governed  by  Maxwell’s 
equations: 

V X H - j + |p  , (1) 


where  the  current  density  j 
time-varying  form  as 


oE  (o  is  the  conductivity).  Equation  (1)  can  be  written  in 


v X H “ (o  + jue) E. 


Since  this  report  deals  with  good  conductors  at  low  frequencies,  cue  is  much  less  than  a 
and  can  be  neglected.  This  assumption  amounts  to  neglecting  displacement  current  in  com- 
parison with  conduction  current.  Thus  one  has 
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and  taking  the  curl  of  each  side  yields 


v X v X H = v(v  • H)  - v2H  ■ 7 X ®E, 


v • B = 0 


H(v  • H)  = 0,  ( 

if  homogeneity  of  permeability  p is  assumed.  Therefore,  substituting  Eq.  (6)  into  Eq.  (4), 
one  has  ’ 

-v2H  = o(v  X E)  ( 


-*■  ■ •(-*)  ■ 


vx*‘~w  ■ 

Hence  the  differential  equation  governing  the  penetration  of  a magnetic  field  into  a con- 
ductor is 


’•■-•Mr. 


where  p is  the  resistivity  and  the  permeability  is  independent  of  time.  Equation  (10)  can  be 
written  in  cylindrical  coordinates  as 


Jaw,  bH.  bH 
p[3ra>-  +"5T  a*  +sraz 

«Ji|5 

p dt  • 


I 


MBBK8 

For  a long  thin  solenoid  with  a rod-shaped  core  and  with  Its  long  axis  in  the  s direction  the 
following  conditions  are  true: 

„ dtf  . a Hr  3H  dHr 

v'"' - S?  - ",  - v2".  - S*  ■ H.  —W  • sr  • I12» 

Thus  Eq.  (11)  reduces  to 

v p s r <«> 

and  with  the  Laplacian  written  in  cylindrical  coordinates  is 


1 ' 1 »*».  8»W,  ,W, 

rvW  7 *r- 


But  for  a long  thin  solenoid 


a2//  a2//, 

1 - o.  ( 

a*2  as2  1 

Thus  the  differential  equation  governing  the  penetration  of  an  alternating  magnetic  field 
into  a rod  contained  within  a long  thin  solenoid  is 


r dr  Vr  c)r  . 


The  solution  of  Eq.  (16)  is  given  by  Scott  [1]  as 


W,(r)  - Ht 


/W\  k 0r\ 


ber  0 + bei  0 


where  r is  the  radial  distance  from  the  cylinder  axis  to  the  point  in  question,  d is  the 
diameter  of  the  rod,  and  Hf  is  the  magnetic  field  intensity  at  the  surface  of  the  rod.  The 
parameter  0 is  given  by 

9 ■ • 1 
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where  p,  p,  and  f are  the  rod’s  incremental  permeability,  resistivity,  and  frequency  respec- 
tively in  the  International  System  of  Units  (SI  units).  The  functions  ber  and  bei  of  Eq.  (17) 
are  the  real  and  imaginary  portions  of  the  zeroth-order  Bessel  function  of  the  first  kind  re- 
spectively and  have  been  tabulated  [2] . 


Define  a characteristic  frequency  fc  by  setting  the  parameter  0 = 1.  Solving  for  f 
yields  c 


(19) 


To  examine  the  physical  significance  otfc,  solve  Eq.  (17)  for  the  magnitude  o(Hz(0)IH 
(that  is,  let  r = 0,  which  is  the  center  of  the  rod)  when  0 = 1,  which  corresponds  to  f - f‘. 
This  yields  Hx(0)/Hg  = 0.985.  Thus  the  magnetic  field  at  the  center  of  the  rod  is  98.5%  of 
its  value  at  the  surface  at’  f - fc  . Therefore  f corresponds  to  the  frequency  for  a given  diam- 
eter,  permeability,  and  resistivity  below  which  the  penetration  of  flux  into  the  rod  is  essen- 
tially uniform. 


The  electrical  impedance  of  the  rod  is 


Z = jojL  = 


ju)N J'  Bt(r ) dA 


(20) 


where  Bg(r)  = pHx(r).  If  one  carries  out  the  integration  indicated  by  Eq.  (20),  the  result  is 


Z = L, 


.If1 


0 bei*  0 — bei  0 ber 
ber2  0 + bei2  0 


’0\  .2/1 


ber  0 ber'  0 + bei  0 bei'  0^ 
ber2  0 + bei2  0 ) 


where 


N^Ap 

~T~ 


, (21) 


(22) 


k 


(24) 
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Eq.  (21)  can  be  rewritten  as 


(25) 


Z - juL0(xn  ~ iXt)  m uL0xf  + ;u>L0x*  . (26) 

The  terms  Xd  and  X/  are  the  real  and  imaginary  parts  of  an  eddy  current  factor  x for  a rod. 

*"  Equation  (26)  can  be  alternatively  written  as  .... 

Z - juL0x  - juLqXqV*,  (27) 


where 


Xo  “ (x*  + 


and 


r ■ 

Thus  the  electrical  impedance  can  be  viewed  as  an  inductor  of  value  LQx.  The  complex  eddy- 
current  factor  x may  be  attached  to  the  permeability  in  LQ,  thus  accounting  for  eddy  currents 
by  viewing  the  permeability  as  complex.  The  derivation  of  xR  and  X/  is  originally  due  to 
Scott  (1] . The  expressions  for  xR  and  X;  can  be  found  in  both  Bozorth  (3]  and  Scott  (1) . 

The  permeability  p'  in  Eq.  (24)  is  the  "apparent  permeability”;  that  is,  it  is  the  perme- 
ability which  the  magnetic  field  sees  in  the  presence  of  eddy  currents.  The  permeability  p in 
Eq.  (24)  is  the  limiting  value  of  p\  approached  at  low  frequencies.  The  A/?  in  Eq.  (25)  is  the 
difference  between  the  measured  alternating-current  resistive  component  of  the  impedance 
R and  the  direct-current  resistance  Rq.  The  difference  A R is  caused  by  eddy  currents.  The 
inductance  L0  in  Eq.  (22)  or  (25)  is  the  limiting  value  of  the  measured  inductance  approached 
at  low  frequencies.  Figure  2 shows  plots  of  Xr  , X/.  Xo<  ***1  ?•  Several  useful  approximations 
to  X* . X/.  »nd  t for  high  and  low  f/fe  are  also  shown  in  Fig.  2.  The  use  of  Eq.  (24)  and  Fig.  2 
shows  that  the  “apparent  permeability”  p'  is  a decreasing  function  of  frequency. 


Fig  2 — Eddy  current  component*  for  a rod 


Prom  Eq.  (26)  it  is  seen  that  in  the  presence  of  eddy  currents  the  electrical  impedance 

Stakes  on  a resistive  component  which  is  proportional  to  The  physical  significance  of  the 
resistive  component  is  that  it  is  the  energy -dissipating  mechanism  of  the  eddy  currents.  It  is 
also  seen  that  the  imaginary  portion  of  the  impedance  is  proportional  to  \„.  This  results  in 
a reduction  in  the  inductance  as  a function  of  frequency.  This  observation  is  consistent  with 
the  decrease  in  the  "apparent  permeability”  as  frequency  increases.  The  equivalent  circuit 
of  the  electrical  impedance  of  a long  thin  solenoid  with  a conducting  rod-shaped  core  in 
the  presence  of  eddy  currents  is  shown  in  the  general  case  in  Pig.  3. 

The  exact  equivalent  circuit  shown  in  Pig.  3 is  inconvenient  to  use,  because  both  K and 
L depend  on  frequency  A much  more  convenient  expression  is  a lumped  equivalent  circuit 
in  which  none  of  the  components  vary  with  frequency.  This  can  be  derived  by  considering 
the  expressions  for  \K  and  given  by  Boaorth  |3)  for  a rod  when  f ^ f \ 
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Pig.  3 — The  exact  equivalent 
circuit  of  a long  thin  aolenoid 
with  a conducting  rod-ehaped 
core  in  the  preaence  of  eddy 
current* 


Substitute  the  fint  term  of  xA  and  X/  into  Eq.  (26).  The  result  is 


Z * 


1 / + >wLo 

# C 


(30) 


for  f/fe  < 1.  This  representation  is  still  somewhat  inconvenient,  since  the  real  term  depends 
on  frequency  squared.  If  the  series  representation  of  Eq.  (30)  is  converted  to  its  equivalent 
parallel  form,*  the  parallel  resistance  and  Inductance  are'  * 


e ••  • • • 


Rp  - 16a/, 


Hi * iS 


(31) 


and 


(32) 


When  f < fe,  the  term  /2/64/2  can  be  neglected  with  an  error  of  less  than  2%.  Thus  the 
lumped  equivalent  circuit  of  a conducting  rod  in  the  presence  of  eddy  currents  shown  in 
Fig.  4 is  obtained.  The  characteristic  frequency  fe  in  Fig.  4 is  given  by  Eq.  (19). 

A derivation  analogous  to  the  one  that  leads  to  Fig.  3 can  be  performed  for  a conduct- 
ing core  composed  of  thin  laminated  sheets.  Only  the  results  will  be  presented  here.  The 
original  derivation  can  be  found  in  Ref.  1.  The  eddy  current  components  for  a thin  sheet 
are  given  by  Scott  [1]  and  Bozorth  [3]  and  are 


(33) 


8 


•‘i  llT*  lJ*~iiiV*  " v 
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Fig.  4 — The  approximate 
equivalent  circuit  of  a 
aolenoid  with  a conducting 
rod-ahaped  core  in  the 
presence  of  eddy  currents 
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In  Eqs.  (33)  and  (34) 


f — 2£- 
c nt  2p 


where  p is  the  resistivity  in  ohm-meters,  t is  the  thickness  in  meters,  and  p is  the  permeabil- 
ity in  SI  units.  The  characteristic  frequency  for  the  case  of  a thin  sheet  again  means  the 
frequency  below  which  the  penetration  of  flux  into  a sheet  is  essentially  uniform.  Plots  of 
Xjft  X/»  X0>  and  f versus  f/fc  can  be  found  in  Ref.  4.  The  exact  equivalent  circuit  of  the  elec- 
trical impedance  of  a long  thin  solenoid  with  a conducting  core  composed  of  thin  lamina- 
tions in  the  presence  of  eddy  currents  is  also  given  by  Fig.  3,  with  Eqs.  (33)  and  (34)  now 
applying  to  xR  and  \r 

An  approximation  to  the  exact  equivalent  circuit  of  a long  thin  solenoid  with  a core 
composed  of  thin  sheets  has  been  obtained  in  Ref.  4.  The  approximate  equivalent  circuit  is 
similar  to  Fig.  4 except  that  the  resistance  is  67 rfcLQ,  where  f is  given  by  Eq.  (35). 

The  complex  impedance  of  the  equivalent  circuit  shown  in  Fig.  4 is 


Zc  =“i— 

bj  + 


167rfoj2Lr 


(16*4) 


- + j — i — 
2 ^ * 


(I6.4)2  IOLq 

2 * (16»4)2 
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Let  the  real  and  imaginary  portions  of  Eq,  (36)  equal  Rc  and  A’c  relatively.  Thus 


wcL0 


/ f v + w 

equivalent  V 77/ 

circuit 


WcL0 


(equivalent  \ c 
(circuit 


The  exact  expressions  for  Rc/ucLQ  and  Xf/wfL0  are  given  by  the  real  and  imaginary  por- 
tions respectively  of  Eq.  (26)  normalized  to  wfL0: 


wcL0 


exact 

theory 


fr/ 

- 7T 


WCL0  'c 

exact 

theory 

Figure  6 is  a comparison  of  the  loci  of  points  predicted  by  Eqs.  (37)  and  (38)  (the 
equivalent  circuit  of  Fig.  4)  and  those  predicted  by  Eqs.  (39)  and  (40)  (the  exact  equivalent 
circuit  of  Fig.  3).  The  loci  of  the  impedance  of  the  lumped  equivalent  circuit  describes  an 
arc  of  a circle  (the  lower  curve  in  Fig.  6)  which  agrees  within  1%  with  the  exact  equivalent 
circuit  up  to  f * f . Deviations  begin  al>ove  f - / (particularly  in  the  resistance),  but  agree- 
ment is  still  fair  at  3f  and  usable  up  to  f - \fc . The  loci  of  the  impedance  of  a lossless  in- 
ductor, if  plotted  on  Fig.  6,  would  appear  as  a straight  line  on  the  ordinate  with  the 
frequency  scale  increasing  upward. 

The  comparison  of  the  exact  and  approximate  equivalent  circuits  for  a core  composed 
of  a stack  of  thin  laminated  sheets  is  given  in  Ref.  4.  In  this  case  the  approximate  equivalent 
circuit  is  usable  up  to  f • 2 ff. 
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Fig.  5 — The  normalized  theoretical  components  of  the  impedance 
of  a core  composed  of  a rod  (upper  curve)  and  the  normalized  com- 
ponenta  of  the  impedance  of  the  equivalent  circuit  shown  (lower 
curve) 


An  enlightening  alternate  derivation  of  the  equivalent  circuit  of  Fig.  4 can  be  obtained 
through  use  of  Faraday’s  law.  Consider  Fig.  1,  which  shows  a conducting  rod  in  the  presence 
of  an  alternating  magnetic  flux.  Assume  that  the  flux  density  is  increasing  and  uniform 
across  the  cross-sectional  area  of  the  rod.  This  assumption  means  that  the  rod  is  very  long 
and  that  the  frequency  is  below  the  characteristic  frequency.  Faraday’s  law  relates  the  time- 
changing magnetic  flux  to  the  electric  field  produced  for  a fixed  region  in  space: 


a ds 


- ( 
dtJA 


B * da. 


11 


f'j 


^ ... . 


(41) 
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Now  perform  the  line  integral  of  the  electric  field  about  the  circular  path  a diitance  r from 
the  axis.  The  electric  field  is  constant  in  magnitude  and  tangential  around  this  path  due  to 
the  circular  symmetry.  Hie  angle  between  E and  ds  is  0°.  The  flux  density  B is  also  uniform 
across  the  rod,  and  the  angle  between  B and  dA  is  0°.  If  one  considers  only  the  magnitude, 
then 


Thus 

E2irr  ~ It  (*r2B)  (43) 


■ wr 


2 dfi 
dt  ' 


(44) 


which  implies  that 


E .1™ 
t 2 dt 


(46) 


The  density  of  the  circumferential  surface  current  at  a radius  r is 


J ml  m JLM 

p 2p  dt 


The  total  instantaneous  power  lost  is 

p (t)  ~^)pJ2  dv 


« 2»  2 


' f f f ‘‘'ffirf '***• 


m-0  t-0  r-0 


which  yields 


( 

\ 


P(0 


wdH  <fa  /dB\* 
4 32 p\dtJ' 
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which  implies  that 

R . J2pA^l  ohms.  (54] 

d2C 

The  resistance  given  by  Bq.  (54)  is  identical  to  the  resistance  I6n ffLQ  shown  in  Fig.  4.  Thus 
both  derivations  yield  the  same  result  for  a long  rod  at  frequencies  below  the  characteristic 
frequency.  The  resistance  GitfeLQ  of  a thin  laminated  sheet  can  be  obtained  in  an  analogous 
manner  through  use  of  Faraday’s  law  [6] . 


v- 


Equivalent  Circuit  of  a Magnetostrictive  Toroid  in  the 
Presence  of  Eddy  Currents 

The  equivalent  circuits  presented  thus  far  have  had  no  mention  of  magnetostriction. 
This  subsection  will  present  the  equivalent  circuit  of  a magnetostrictive  toroid  with  a core 
composed  of  a circularly  cross-sectioned  rod  or  a number  of  thin  laminated  sheets  in  the 
presence  of  eddy  currents. 

The  equivalent  circuit  of  a magnetostrictive  transducer  can  be  formed  in  three  ways: 
• by  the  mobility  analogy, 


13 


MEEKS 


• by  shifting  force  and  velocity  by  90°,  or 

• by  F.V.  Hunt’s  space  operator. 

The  effect  of  eddy  currents  on  the  equivalent  circuit  formed  by  shifting  the  force  and  veloc- 
ity by  90°  is  discussed  in  Refs.  4 and  6.  The  effect  of  eddy  currents  on  the  equivalent 
circuit  formed  by  Hunt’s  space  operator  is  discussed  in  Ref.  7 along  with  an  enlightening 
discussion  of  all  three  representations  in  the  absence  of  eddy  currents.  The  equivalent 
circuits  formed  by  shifting  the  force  and  velocity  by  90°  and  formed  by  Hunt’s  space 
operator  have  the  disadvantage  of  an  imaginary  frequency-dependent  turns  ratio,  which 
leads  to  frequency -dependent  mechanical  elements  even  in  the  absence  of  eddy  currents. 

The  presence  of  eddy  currents  further  complicates  these  two  equivalent  circuits.  The  mobil- 
ity analogy  is  the  simplest  to  use  and  the  most  widely  accepted  representation  of  electro- 
magnetic transducers.  The  use  of  the  mobility  analogy  when  eddy  currents  are  present 
leads  to  an  equivalent  circuit  that  has  frequency-independent  elements  for  a certain  fre- 
quency range,  as  will  be  shown. 

The  equivalent  circuit  of  a magnetostrictive  toroid  in  the  mobility  analogy  without 
eddy  current  loss  is  shown  in  Fig.  6.  This  equivalent  circuit  does  not  include  ohmic  loss, 
which  can  be  accounted  for  by  simply  inserting  a resistor,  whose  value  is  equal  to  the  dc 
resistance  of  the  coil,  in  series  with  the  blocked  inductance  L0.  In  the  mobility  analogy, 
force  is  analogous  to  current,  and  velocity  is  analogous  to  voltage.  In  this  figure,  /,0  is  the 
blocked  or  purely  electrical  impedance,  Cm  is  the  mechanical  compliance,  A/m  is  the  me- 
chanical mass,  1/Rm  is  the  inverse  of  the  mechanical  resistance,  and  1 /RA  and  MA  compose 
the  radiation  mobility  of  a sphere.  The  toroid  in  this  report  is  assumed  to  be  end  capped, 
so  that  it  is  a volume  expander.  At  low  frequencies  the  radiation  mobility  of  a sphere  is 
used  to  approximate  the  toroid’s  radiation  mobility.  The  expressions  for  the  equivalent 
circuit  parameters  are 


L 


0 


N2vf3 

2irr 


(55) 


Force  Pressure 


Fig.  6 — The  equivalent  circuit  of  a magnetostrictive  toroid  in  the 
mobility  analogy  without  eddy  current  loaa 
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(60) 


where 

N is  the  number  of  turns  of  wire, 

Ac  is  the  cross-sectional  area  of  the  metallic  core, 

/if  3 is  the  dynamic  permeability  of  the  core  at  a constant  strain  (the  blocked  perme- 
ability), 

r is  the  radius  of  the  toroid, 

D is  the  density  of  the  core, 

s"  is  the  elastic  compliance  coefficient  measured  at  constant  magnetic-field  intensity, 
a>0  is  2jt  times  the  unloaded  (in  air)  resonance  frequency, 

Qm  is  the  mechanical  Q measured  in  air, 

rf  is  the  equivalent  radius  of  a sphere  which  has  a surface  area  equal  to  the  active  sur- 
face area  of  the  end-capped  toroid, 

dQ  is  the  density  of  water, 

cQ  is  the  speed  of  sound  in  water, 

0 is  the  electromechanical  turns  ratio,  and 
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A is  the  mechanoacoustic  turns  ratio,  which  in  this  case  is  the  active  surface  area  of 
the  toroid. 

Eddy  currents  are  accounted  for  in  the  equivalent  circuit  by  assuming  a complex 
permeability  as  suggested  by  the  LQx  term  of  Eq.  (27).  The  eddy  current  factor  x appears 
not  only  in  the  blocked  impedance  but  also  in  the  electromechanical  turns  ratio  <f>  in  the 
mobility  analogy,  as  will  be  shown. 

The  derivation  of  the  electromechanical  turns  ratio  0 begins  with  the  magnetostrictive 
equations  of  state  for  a thin  rod: 


S3  " *33^3  + *33B3 
and 

H3m~*33T3+v33B3  • 


(61 


(62 


If  S3  is  set  equal  to  zero,  Eq.  (61)  becomes 


T 

3 s®  ' 
*33 

The  flux  density  in  the  clamped  (S3  ■ 0)  toroid  is 

B3  “ ^33x^3- 


(63 


(64 


Here  the  permeability  is  complex  due  to  the  presence  of  eddy  currents.  Substitution  of 
Eq.  (64)  into  Eq.  (63)  yields 


~*33*43X^3 

s® 

*33 


But  Ts  - ; therefore 


F3  “ 


^33^33  *AeH3 

gB 

*33 


(66 


(66 


In  a toroid  the  radial  force  is  related  to  the  circumferential  force  in  the  toroid  by 


the  following  equation  is  produced  for  magnitude  only: 


£r  Nd33Ac* 

/ * rJH 
rs33 


0 


X • 


without 
eddy  currents 


(73) 


Since  the  eddy  current  factor  x is  a complex  number,  the  turns  ratio  Ff/I  in  Eq.  (73)  is  also 
a complex  number. 


The  equivalent  circuit  of  Fig.  6 can  now  be  modified  to  account  for  eddy  current  losses 
by  replacing  the  blocked  inductance  LQ  by  LQx  and  the  electromechanical  turns  ratio  0 by 
0X-  Figure  7 shows  the  equivalent  circuit  of  a magnetostrictive  toroid  in  the  presence  of 
eddy  currents  with  the  transformers  removed.  The  complex  eddy-current  tactor  is  reflected 
in  the  mechanical  mobility  and  the  radiation  mobility.  If  the  individual  components  of 
Fig.  7 are  separated  into  their  real  and  imaginary  portions,  Fig.  8 results.  Each  of  the  com- 
ponents in  Fig.  8 is  a complicated  function  of  frequency , due  to  the  eddy  current  terms.  This 
makes  the  frequency  analysis  of  the  circuit  difficult. 

To  make  the  analysis  of  the  circuit  somewhat  simpler,  it  was  necessary  to  convert  the 
circuit  to  mechanical  units.  This  is  done  by  moving  the  electromechanical  transformer  of 
Fig.  6 to  the  left  and  the  mechanoacoustic  transformer  to  the  right.  The  result  is  shown  in 
Fig.  9.  The  resistors  in  Fig.  9 represent  inverse  mechanical  resistances  with  units  of  seconds 
per  kilogram.  The  inductors  represent  mechanical  compliances  with  units  of  meters  per 
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Flf.  8 — Hi*  equivalent  circuit,  in  electrical  unit*,  of  a magneto*  trictive  toroid  in  the 
peace  nee  of  eddy  current*,  with  individual  element*  separated  into  real  and  imaginary 
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Fig.  9 — The  equivalent  circuit,  in  mechanical  units,  of  a magnetostrictiive 
toroid  in  the  presence  of  eddy  currents 


newton.  The  capacitors  represent  masses  with  unite  of  kilograms.  The  through  variable  is 
the  force  produced  by  the  toroid  in  newtons.  The  across  variable  is  the  velocity  of  the 
toroid  in  meters  per  second.  The  eddy  current  terms  in  Pig.  9 are  lumped  into  three  com- 
ponents: the  two  elements  of  the  blocked  mechanical  mobility  (branches  2 and  3 (B2  and 
B3))  and  the  force  generator  (branch  1).  The  ideal  sinusoidal  current  generator  of  Pig.  8 
has  now  become  a force  generator  of  value  I<f>x  newtons.  Figure  2 shows  that  x0  and  thus 
the  magnitude  of  the  force  generator  are  decreasing  functions  of  frequency.  The  decrease  in 
force  produced  by  the  magnetostrictive  core  as  a function  of  frequency  is  due  to  the  reverse 
flux  generated  by  the  eddy  currents.  The  monotonic  decrease  of  force  with  frequency  may 
be  explained  by  noting  that  the  eddy  currents  increase  with  frequency.  Hence  the  opposing 
magnetic  field  increases,  thereby  resulting  in  a continual  decrease  in  the  net  flux  that  pro- 
duces the  force  from  the  magnetostrictive  core. 

The  equivalent  circuit  shown  in  Fig.  9 can  now  be  used  to  analyze  a magnetostrictive 
transducer  in  the  presence  of  eddy  currents.  Figure  9 can  be  used  to  calculate  the  efficiency, 
transmitting  responses,  receiving  response,  and  impedance,  all  as  a function  of  frequency. 
Before  performing  the  analysis,  one  must  fiwt  observe  how  the  power  dissipated  in  the 
circuit  elements  of  Fig.  9 is  calculated. 

Power  is  defined  as 


P-Re(/*e),  (74) 

where  /*  is  the  complex  conjugate  of  the  current  and  e is  the  voltage.  Figure  9 implies  that 


I* 


(75) 
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where  V is  the  velocity  and  F is  the  force.  Thus 


therefore 


where  Ym  is  the  mechanical  admittance  (mobility)  whi 
mechanical  impednace.  Equation  (79)  can  be  written  as 


for  a rod  of  circular  cross  section  and  by 


symbols  in  Eqs.  (81)  and  (82)  are  the  same  as  in  Eqs.  (19)  and  (35) 
i33  is  the  blocked  incremental  permeability.  The  blocked  characteristic 
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frequency  is  required  whenever  one  is  dealing  with  an  equivalent  circuit,  since  the  eddy  cur- 
rent factor  (which  is  written  as  a series  composed  of  f/f*  as  in  Eqs.  (28)  and  (29))  is 
attached  to  the  blocked  permeability,  as  in  Eq.  (64),  and  to  the  blocked  inductance  of 
Fig.  7.  Once  f/f*  is  determined,  Xq  is  read  from  Fig.  2 or  from  Ref.  4 if  a thin  sheet  is  being 
used.  The  magnitude  of  the  force  generator  is  thus  IQXq-  The  phase  angle  of  x is  not  put  into 
the  force  generator,  since  one  is  interested  in  only  the  magnitude  of  the  velocity,  as  indi- 
cated in  Eq.  (80).  The  analysis  of  Fig.  9 was  carried  out  on  a PDP-11/45  computer  using  an 
ac -circuit  analysis  program  (8J . The  values  of  the  elements  were  chosen  to  be  numerically 
equal  to  their  mechanical  quantities.  The  branch  voltages  and  currents  were  therefore  inter- 
preted as  velocities  and  forces  respectively.  The  power  dissipated  in  each  element  was  cal- 
culated using  the  velocity  from  the  circuit  analysis  program,  the  phase  angle  f from  Fig.  2 
or  Ref.  4,  and  Eq.  (80).  It  is  not  necessary  to  calculate  the  power  dissipated  in  the  blocked 
mechanical  mobility  (branches  2 and  3 of  Fig.  9)  using  the  method  described.  The  power 
dissipated  in  these  elements  is 

P • l2uLQxr  (83) 

which  shows  that  this  power  can  be  determined  from  known  constants  and  a given 
frequency. 

The  steps  necessary  to  perform  an  efficiency  analysis  using  Fig.  9 are  summarized  as 
follows: 

1.  Calculate  f/f*  (to m the  desired  frequency  and  Eqs.  (81)  or  (82). 

2.  Determine  x^  from  Fig.  2 or  Ref.  4. 

3.  Input  the  value  of  the  force  generator  (/0Xq)  and  the  values  of  mechanical  com- 
pliances, masses,  and  resistances  into  an  ac-circuit  analysis  program. 

4.  Obtain  the  magnitude  of  the  velocity  across  each  element  (except  the  blocked 
mechanical  mobility)  and  the  phase  angle  ? from  the  graphs  to  calculate  the  power 
dissipated  in  each  element  via  Eq.  (80). 

5.  Use  the  power  dissipated  in  the  mechanical  elements  and  Eq.  (83)  to  determine  the 
efficiency  in  the  presence  of  eddy  currents: 

PA 

efficiency  - 4 Pr  Vp!  • (84) 

M A A 

where  PA  is  the  acoustic  power,  Py  is  the  mechanical  power,  and  PE  is  the  power 
dissipated  by  eddy  currents. 

6.  To  calculate  the  efficiency  at  another  frequency,  return  to  step  1 and  use  the  new 
f/f*  to  calculate  the  new  magnitude  of  the  force  generator.  This  new  value  of  the 
force  generator  is  input  into  the  ac-circuit  analysis  program,  and  the  new  power 
dissipated  in  each  element  is  calculated  from  Eq.  (80).  For  each  new  frequency  only 
one  element  (the  force  generator)  of  Fig.  9 need  be  changed  to  calculate  the 
efficiency. 
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To  calculate  the  transmitting  current  response,  one  first  calculates  the  acoustic  power 
dissipated  in  the  mechanical  radiation  mobility  (1  /A  2RA  ) using  the  technique  described. 

*“  ▲! * a.  ? • * > . > . 


Then  the  on-axis  acoustic  pressure  is  calculated  via 

„ . i&W.t 

P r \ 4ir 


where  p is  the  on-axis  acoustic  pressure,  Pq  is  the  acoustic  power,  d»  is  the  density  of  water, 
c0  is  the  speed  of  sound  in  water,  and  Rg  is  the  directivity  factor.  Since  the  pressure  cal- 
culated is  for  a 1-ampere  drive  (if  the  I in  /0Xo  is  set  equal  to  unity),  the  pressure  obtained 
from  Eq.  (85)  will  be  the  transmitting  current  response.  The  free-field  voltage  sensitivity  can 
be  obtained  [9]  from 

M 

J = reciprocity  parameter  (86) 

where  M is  the  free-field  voltage  sensitivity  and  S is  the  transmitting  current  response, 
along  with 

•/=VX1°12’  (87) 

in  which  d is  the  measurement  distance,  dQ  is  the  density  of  water,  and  f is  the  frequency, 
all  in  SI  units. 

To  calculate  the  electrical  impedance,  one  proceeds  as  follows: 

V e 1 ze 


e 1 

Z' 

(88) 

0X/0X 

02x2 

V ,2  2 
= 70V 

V . 

= -7-0X, 

(89) 

where  V is  the  velocity  and  / is  the  current.  Equation  (89)  can  be  simplified  for  calculation 
by  noting  that  the  velocity  will  contain  the  phase  angle  — f when  the  force  generator  has  the 
angle  — If  one  uses  only  the  magnitude  of  the  force  generator  in  the  equivalent  circuit  and 
adds  — f to  the  phase  angle  of  the  resulting  velocity,  then 


-70X0  e~J2t 


and  if  / = 1 ampere  and  V = I V I e*9  , then 

Zt  * IVl0XoC><9-2f)  . ( 

The  procedure  for  determining  Zf  from  the  equivalent  circuit  of  Fig.  9 can  be  sum- 
marized as  follows: 

1.  Follow  steps  1,  2,  and  3 that  were  listed  for  the  efficiency  analysis. 
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2.  Calculate  the  values  of  the  blocked  mechanical  mobility  for  the  particular  frequency 
and  ratio  of  f to  f*.  Input  these  values  into  the  ac -circuit  analysis  program. 


3.  Record  the  magnitude  I Viand  phase  0 of  the  input  (node  1,  Fig.  9)  velocity.  Use 
Eq.  (91)  to  calculate  Zf. 

4.  To  calculate  Zf  at  another  frequency,  return  to  step  1.  Notice  that  three  circuit 
elements  must  be  changed  for  each  frequency  point.  These  are  the  two  blocked 
mechanical  mobilities  and  the  force  generator. 

The  transmitting  voltage  response  is  obtained  as  follows: 

1.  Follow  steps  1 and  2 of  the  procedure  just  summarized  for  calculating  Zg. 

2.  Record  the  magnitude  I V I of  the  input  velocity  (node  1,  Fig.  9),  and  use  the  follow- 
ing equation  to  calculate  the  input  voltage: 

lei  = magnitude  of  the  voltage  = IVI^Xq-  (92) 

3.  Calculate  the  power  dissipated  in  the  radiation  mobility  (1  /A2Ra  ) for  the  input 
voltage  of  step  2.  Use  Eq.  (85)  to  calculate  the  on-axis  acoustic  pressure  and  divide 
this  pressure  by  the  result  of  Eq.  (92). 

The  equivalent  circuit  of  Fig.  9 is  based  on  an  exact  eddy-current  theory  that  will 
work  at  any  frequency  subject  to  the  assumptions  under  which  the  circuit  was  derived.  But 
the  wide  frequency  range  of  the  equivalent  circuit  of  Fig.  9 is  a tradeoff  with  the  com- 
plexity of  the  analysis  required  to  obtain  results.  It  would  be  desirable  to  have  an  equivalent 
circuit  that  has  lumpe  ' components,  that  is,  with  none  of  the  elements  (capacitors,  resistors, 
inductors,  or  current  generator)  changing  with  frequency.  Such  an  equivalent  circuit  is  pos- 
sible over  a limited  frequency  range  and  will  be  described  in  the  succeeding  paragraphs. 

To  derive  a lumped  equivalent  circuit,  one  must  approximate  the  eddy  current  factor  x 
in  such  a way  that  the  frequency-dependent  elements  of  Fig.  7 become  a combination  of 
frequency-independent  elements.  This  has  already  been  accomplished  for  the  blocked  im- 
pedance (branch  2)  of  Fig.  7.  When  f/f*  < 2,  the  blocked  impedance  (for  a circularly  cross- 
sectioned  core)  of  Fig.  7 can  be  represented  quite  accurately  by  Fig.  4,  where  LQ  is  given  by 
Eq.  (55)  and  f*  is  given  by  Eq.  (81).  Figure  4 was  derived  by  assuming  xR  ^ 1 and 
Xt  58  f/8f*,  which  implies  that 


X = XR  ~ JXi 


(93) 


(94) 
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for  f/fl  < 2.  Therefore  branch  3 of  Fig.  7 become* 


juCm<t>2x2  » iucm<t>2 


- 2/ — 


1 - 2/- 2 

fy;  64 


* j“Cm<t>2  *~2jf  — • 

if  the  last  term  of  Eq.  (96)  is  neglected.  Equation  (96)  can  be  written  in  an  equivalent 
parallel  representation  shown  in  Fig.  10.  Figure  10  is  the  lumped  form  which  is  desired. 
Branch  4 in  Fig.  7 becomes 


(95) 
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Equation  (98)  can  be  represented  by  Fig.  11.  Branch  6 in  Fig.  7 becomes 
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Fig.  10  — Approximate 
equivalent  circuit  of  the 
mechanical  compliance 
in  the  pretence  of  small 
eddy  current* 
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Fig.  11  — Approximate  equivalent  circuit 
of  the  mechanical  maaa  in  the  ^reaence  of 
•mall  eddy  current! 


Equation  (100)  can  be  represented  by  Fig.  12.  Branches  6 and  7 are  represented  similarly 
The  resulting  lumped  equivalent  circuit  of  a magnetostrictive  toroid  (with  a circular-cross- 
section core)  in  the  presence  of  eddy  currents  (f/f*  < 2)  is  shown  in  Fig.  13.  Figure  13  has 
nine  nodes  and  15  branches  compared  with  three  nodes  and  seven  branches  for  the  equiva- 
lent circuit  without  eddy  currents.  Figure  13  is  in  electrical  units.  The  efficiency  of  a trans- 
ducer described  by  Fig.  13  is  obtained  by  analyzing  the  circuit  with  an  ac -network  routine 
The  power  dissipated  in  the  resistances  of  Fig.  13  is 


(101) 


where  e is  the  branch  voltage  and  R is  the  branch  resistance.  The  power  dissipated  by  eddy 
currents  is  that  dissipated  in  branches  3,  4,  7,  and  14.  The  power  dissipated  in  the  negative 
resistances  of  branches  7 and  14  must  be  interpreted  as  negative  power  in  order  to  obtain 
correct  results.  The  mechanical  dissipation  occurs  in  branch  9,  and  the  acoustical  dissipa- 
tion  occurs  in  branch  11.  The  efficiency  is  given  by  Eq.  (84).  As  eddy  currents  become 
small,  approaches  infinity,  which  means  that  the  elements  of  Fig.  13  approach  those  of 
the  equivalent  circuit  without  eddy  currents. 

The  transmitting  current  response  is  calculated  from  Eq.  (85)  using  the  power  dis- 
sipated in  branch  11  for  a 1-ampere  current  generator.  The  electrical  impedance  is  calculated 
from  the  quotient  of  the  complex  voltage  at  node  1 and  the  value  of  the  current  generator. 
The  transmitting  voltage  response  is  computed  from  the  power  dissipated  in  branch  11  and 
the  driving  voltage  at  node  1 . 

An  equivalent  circuit  similar  to  Fig.  13  can  be  derived  for  a magnetostrictive  toroid 
with  a core  composed  of  thin  sheets  in  the  presence  of  small  eddy  currents.  In  the  case  of 
a thin-sheet  core  the  eddy  current  factor  x becomes 

X = X«  ~ iX,  (102) 


drcuTt^^T  ^ Eq'  (82)‘  APP‘ying  Eq- (103)  t0  Fig-  7 yie,ds  the  de8ired  equivalent 
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Results  and  Discussion  of  Equivalent  Circuits 

Figure  15  shows  the  electroacoustic  efficiency  in  decibels  versus  normalized  frequency 
and  f/f*  for  a toroid  with  a circular  cross  section.  This  figure  was  generated  from  Figs.  13 
and  9 lor  the  small-eddy-current  approximation  and  the  exact  eddy-current  theory  respec- 
tively. Below  f/f*  - 1 the  exact  and  approximate  theory  are  in  excellent  agreement.  The 
difference  between  the  exact  and  approximate  theory  does  not  reach  1 dB  until  f/f*  “ 2.3. 
The  point  of  maximum  efficiency  need  not  occur  at  the  point  of  maximum  power  output 
f/fr  » 1,  since  here  maximum  efficiency  occurs  at  f/fr  m 1.3. 

Figure  16  shows  the  electroacoustic  efficiency  in  decibels  versus  normalized  frequency 
and  f/f*  for  a toroid  with  a core  composed  of  thin  sheets.  This  graph  was  generated  from 
Figs.  14  and  9 for  the  small-eddy -current  approximation  and  the  exact  eddy-current  theory 
respectively.  The  agreement  between  the  exact  and  approximate  theories  is  not  as  good  here 
as  in  Fig.  16.  There  is  a 1-dB  difference  at  f/f*  ■ 0.75  and  1.7-dB  difference  at  f/f*  ■ 1.0. 
Therefore  the  use  of  Fig.  14  for  efficiency  calculations  should  be  limited  to  ratios  of  f to 
f*  below  unity.  The  predictions  of  Figs.  13  and  14  yield  a greater  value  of  efficiency  than 
would  be  measured. 


Fl«.  16  — The  electroacoustic  efficiency  of  a magnetoetrictive  toroid 
with  a circular  crow  section 
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Fig.  16  — The  electroacoustic  efficiency  of  a magnetoatrictive  toroid 
with  a core  composed  of  thin  sheets 


Figure  17  shows  the  electrical  impedance  versus  frequency  for  a magnet .ostrictive 
toroid  with  a circular  cross  section.  This  figure  was  generated  from  Figs.  13  and  9 for  the 
small -eddy-current  approximation  and  the  exact  eddy-current  theory  respectively.  The 
approximate  and  exact  theories  show  good  agreement  over  the  frequency  range  covered. 
The  approximate  equivalent  circuit  of  Fig.  13  should  give  useful  impedance  information 
up  to  f/f*  “ 2.  Figure  14  can  be  used  to  yield  similar  results  for  a thin-sheet  core,  although 
it  is  limited  to  f/f*  < 1.  A practical,  highly  efficient  magnetostrictive  transducer  is  designed 
so  that  f/f*  < 2 over  its  operating  bandwidth.  Therefore  the  approximate  equivalent  circuits 
of  Figs.  13  and  14  are  quite  useful  from  a practical  point  of  view. 

Figure  18  shows  the  transmitting  current  response  for  a theoretical  magnetostrictive 
toroid  with  a circular-cross-section  core.  This  figure  compares  the  small-eddy-current 
approximation,  the  exact  eddy-current  theory,  and  the  equivalent  circuit  without  eddy  cur- 
rents. The  equivalent  circuit  without  eddy  currents  and  the  small-eddy-current  approxima- 
tion show  little  difference  over  the  frequency  range  considered.  Below  f/f*  - 2 there  is  less 
than  1-dB  difference  between  the  no-eddy-current  results  and  the  exact  eody-current 
theory.  Therefore,  to  calculate  the  transmitting  current  response  at  frequency  ratios  below 
f/f*  - 2,  one  need  use  only  the  equivalent  circuit  without  eddy  current  loss  or  the  small- 
eddy-current  approximation,  since  there  is  little  difference  between  them.  Above  f/f*  - 2 
one  should  use  the  exact  eddy-current  theory.  f 
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Fig.  18  — Transmitting  current  rmponw  (TCR)  for  a magnetoatrictive  toroid 
with  a circular  Croat  arc t ion 


Figure  19  shows  the  free-field  voltage  sensitivity  for  the  same  magnetostrictive  trans- 
ducer as  Fiv.  18.  Again  this  figure  shows  that  below  f/f*  - 2 one  need  use  only  the  equiva- 
lent circuit  without  eddy  currents  or  the  small-eddy-current  approximation  to  calculate 
the  FFVS.  Above  f/f*  - 2 one  needs  to  use  the  exact  eddy-current  theory  to  obtain  good 
results. 

Figure  20  shows  the  transmitting  voltage  response  versus  frequency  for  the  same  mag- 
netostrictive toroid  as  Figs.  18  and  19.  The  equivalent  circuit  without  eddy  current  loss 
(Fig.  6)  is  only  1 dB  above  the  exact  eddy-current  theory  at  f/f*  • 3.  Thus  f/f* » 3 is  the 
dividing  line  between  using  the  exact  eddy-current  theory  and  the  equivalent  circuit  without 
eddy  current  loss.  In  this  case  the  small-eddy-current  approximation  is  not  quite  as  good  as 
the  results  obtained  from  the  equivalent  circuit  without  eddy  currents. 
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Fig.  19  — Free-field  voltage  aenaitivity  (FFV8)  for  the  tranaducer  of  Fig.  18 
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Fig.  20  — Traiumitting  voltage  reapoiue  (TVR)  for  the  traruducer  of  Fig.  18 


The  Effect  of  Magnetomechanical  Coupling  on  Eddy  Current  Loss 

Figure  21  shows  the  equivalent  circuit,  in  electrical  units,  of  an  unloaded  magnetostric- 
tive  toroid  in  the  presence  of  eddy  currents.  This  section  will  show  how  magnetoelastic 
effects  alter  the  amount  of  power  dissipated  by  eddy  currents. 

If  one  calculates  the  power  dissipated  by  eddy  currents  in  branches  2,  3,  and  4 (no 
eddy  current  power  being  dissipated  in  branch  5,  only  mechanical  power),  the  following 
equation  results: 


P.ddy<w)  " X/ 


1 +. 


(104) 
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Pig.  21  — Hi*  equivalent  circuit,  in  electrical  unite,  of 
an  unloaded  magnetoetrictive  toroid  in  the  pretence  of 
eddy  currents 


where 


The  frequencies  w,,  and  u> . are  the  resonance  frequency  and  a damping  frequency  respec- 
tively. Equation  (106)  implies  that  as  the  mechanical  damping  Rm  approaches  zero,  the 
damping  frequency  approaches  infinity.  The  result  of  using  the  small-eddy-current  approx 
imations  of  Eq.  (94)  in  Eq.  (104)  is 
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Fig.  22  — The  power  diuipated  by  eddy  current*  a*  a function  of  frequency 

when  ■ 0.45 
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Thus  Eq.  (108)  becomes 


/■*(!  — *33 ) * free  characteristic  frequency. 
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Equation  (111)  is  the  same  as  the  expression  for  the  eddy  current  loss  of  a nonmagnetostric 
tive  toroid  with  permeability  Pg3  for  f < fj.  Equation  (111)  means  that  at  frequencies  tax 
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below  resonance  (where  the  motional  impedance  is  small)  the  power  lost  in  eddy  currents  of 
a magnetostrictive  toroid  is  the  same  as  in  a nonmagnetostrictive  toroid  with  permeability 
H = pj3.  When  co  = u>0,  Eq.  (107)  becomes 


Peddy 


i2u2olo 

16< 


(112) 


Equation  (112)  is  the  same  as  the  expression  for  the  loss  of  a nonmagnetostrictive  toroid 
with  permeability  p = p|3  at  co  = co0  and  / < /* . When  co  » coQ,  Eq.  (107)  becomes 


*eddy  ^ 


jco2Lt(1  - fe|3)2 

16*fT 


72cj2L0 

16^ 


(113) 


which  is  the  same  as  Eq.  (112)  except  that  co  has  replaced  coQ.  Equation  (113)  means  that 
far  above  resonance  the  = 0.45  curve  of  Fig.  22  is  less  than  the  k 33  = 0 curve  by  a factor 

°f (i  - *i3)2- 

The  most  interesting  phenomenon  shown  in  Fig.  22  is  the  10-dB  drop  below  the  k33  =0 
curve  in  the  eddy  current  loss  at  f/f0  = 1.11,  just  past  the  resonance  frequency.  Savage  and 
Abbundi  [10]  have  observed  this  effect  experimentally.  This  present  theory  explains  that 
effect.  Physically  this  drop  in  the  eddy  current  power  dissipated  means  that  the  skin  depth 
of  the  magnetostrictive  material  increases  over  a narrow  frequency  range  just  above  the  reso- 
nance frequency.  It  has  been  found  that  the  depth  of  the  minimum  at  f/f0  - 1.11  is  con- 
trolled by  both  the  damping  frequency  cod  and  the  coupling  coefficient  k33.  The  damping 
frequency  is  related  to  the  mechanical  resistance  Rm  via  Eq.  (106).  A small  mechanical 
damping  (large  wd)  leads  to  a great  drop  in  the  eddy  current  loss  just  past  <o0.  A large 
mechanical  damping  (small  ud)  leads  to  a shallow  drop  in  the  eddy  current  power  above 
resonance.  The  occurance  of  such  a drop  in  eddy  current  loss  could  lead  to  a dramatic  rise 
in  the  efficiency  of  a transducer  that  has  a large  portion  of  its  loss  in  eddy  current  ohmic 
loss. 


Figure  23  shows  the  eddy  current  power  dissipated  as  a function  of  frequency  for 
L0  = 7.36  H,  f*  = 1000  Hz,  ft33  - 0.3,  u>0  = (2ir)100  Hz,  and  cod  = (2w)875  Hz.  The 
reduced  coupling  coefficient  moves  both  the  maximum  and  the  minimum  closer  to  fQ. 

The  maximum  and  minimum  occur  at  0.95  fQ  and  1.055  fQ  respectively.  The  difference 
between  the  curves  above  resonance  is  less  than  that  in  Fig.  22.  This  is  expected  in  light  of 
Eq.  (113). 


Extension  of  the  Small-Eddy -Current  Approximations  to 
Additional  Equivalent  Circuits 

The  only  equivalent  circuit  discussed  so  far  has  been  that  of  a magnetostrictive  toroid. 
The  small-eddy-current  approximations  presented  in  Figs.  13  and  14  apply  to  any  magneto- 
strictive device  with  an  active  core  that  is  circular  in  cross  section  or  composed  of  a stack  of 
thin  sheets.  Figure  24  shows  an  example  in  which  the  approximations  of  Fig.  13  have  been 
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applied  to  the  equivalent  circuit  of  a thin  magnetostrictive  rod  of  circular  cross  section  with 
one  end  clamped.  Branches  1 through  4 are  the  blocked  impedance.  Branches  5 through  9 
are  Mason’s  approximations  to  the  mechanical  admittance  with  the  approximations  of 
Fig.  13  applied.  Branches  10  and  11  are  the  reciprocal  of  the  mechanical  resistance  with  the 
small-eddy-current  approximation  applied.  In  Fig.  24  the  constants  are  given  by 


^0  - 


c_  = 


and 


iy  ^33 

fi 

(114) 

Nd33A 

*33  8 

(115) 

*338 

A ' 

(116) 

= DAS, 

(117) 

w0  Mm 

20 

(118) 

where  6 is  the  length  of  the  rod,  A is  its  cross-sectional  area,  and  D is  its  density,  all  in  SI 
units.  Several  restrictions  apply  to  Fig.  24: 

• One  end  of  the  rod  is  rigidly  clamped,  and  the  other  end  is  free  to  move; 

• The  frequency  is  near  the  quarter- wave  resonance  of  the  rod; 

• The  rod  is  long  and  thin,  so  that  any  demagnetizing  effects  are  small,  or  equivalently 
a well-laminated  low-reluctance  magnetic  return  path  is  provided; 

• The  frequency  range  of  interest  is  less  than  or  equal  to  twice  the  clamped  characteristic 
frequency  (f*). 

SUMMARY  AND  CONCLUSIONS 

The  equivalent  circuit  of  a magnetostrictive  transducer  in  the  presence  of  eddy  currents 
has  been  derived  using  an  exact  eddy-current  theory.  The  exact  eddy-current  theory  can  be 
Hme  consuming  to  use  but  can  be  reduced  to  an  approximate  theory  which  is  simple  and 
'curate  for  frequency  ranges  of  practical  interest.  The  approximate  eddy-current  theory  is 
useful  in  predicting  the  efficiency,  the  impedance,  the  transmitting  responses,  and  the  re- 
ceiving response  in  a range  of  frequencies  near  the  clamped  characteristic  frequency  f*.  At 
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frequencies  greater  than  twice  the  characteristic  frequency  the  exact  eddy-current  theory 
will  yield  accurate  results.  The  approximate  eddy-current  theory  should  have  immediate 
application  by  transducer  designers  working  with  magnetostrictive  transducers. 

It  was  found  that  the  eddy  current  theory  predicted  a sharp  drop  in  the  eddy-current 
loss  in  a narrow  frequency  range  just  past  resonance.  This  effect  was  ex perimen tally 
observed  by  Savage  and  Abbundi  ( 10) . 
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